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High-energy neutrinos may resonate with relic background neutralinos to form short-lived sneu- 
trinos. In some circumstances, the decay chain that leads back to the lightest supersymmetric 
particle would yield few-GeV gamma rays or charged-particle signals. Although resonant coanni- 
hilation would occur at an appreciable rate in our galaxy, the signal in any foreseeable detector is 
unobservably small. 
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The possibility of detecting relic neutrinos by observ- 
ing the resonant annihilation of extremely energetic neu- 
trinos on the background neutrinos through the reaction 
uu -^ Z^ has been the object of extensive studies QiQimI- 
Given the small neutrino masses (< 1 eV) indicated by 
current experimental constraints, suitably intense sources 
of extremely energetic (10^^ - 10^^-eV) cosmic neutrinos 
are required. The positions and shapes of the absorption 
lines in the extremely high-energy neutrinos spectrum 
are influenced by Fermi motion of the relics and by the 
thermal history of the Universe. 

Relic neutrinos have a special standing: According to 
the standard cosmology, neutrinos should be the most 
abundant particles in the Universe, after the photons of 
the cosmic microwave background, provided that they 
are stable over cosmological times. But the weight of 
cosmological observations argues that neutrinos are not 
the only undetected relics. According to the Wilkinson 
Microwave Anisotropy Probe (WMAP) three-year anal- 
ysis 01, the matter fraction of the present Universe is 
n,nh'^ = 0.127+^;[3';'I, where h = 0.73 ±0.03 is the reduced 
Hubble parameter ^] and 17^ is the ratio of the matter 
density to the critical density Qc = 3ff^/87rGN- (Here H 
is the Hubble parameter and Gn is Newton's constant.) 
The baryonic fraction is Qbh^ = 0.0223to;°°og, and the 
neutrino fraction f^^/i^ = (E* "^i^.) /94 eV w 0.0072 < 
Q.bh?' . Accordingly, we have reason to believe that the 
matter fraction is dominated by an unseen "dark-matter" 
component. A popular hypothesis, realized in many ex- 
tensions to the standard model including supersymme- 
try, holds that the dark matter is dominated by a single 
species of weakly interacting massive particle (WIMP). 

There is no confirmed observation of the passage of 
WIMPs through a detector |a, Q . A positive signal re- 
ported by the dark matter experiment DAMA [31 seems 



z/,£ 





.0 
u A 

FIG. 1: Resonant sneutrino formation in neutrino-neutralino 
collisions. Double lines denote superpartners. 



in conflict with upper limits set by the Cryogenic Dark 
Matter Search and the EDELWEISS experiment [Tol |. 
Accelerator experiments have so far not yielded evidence 
for the production of a superpartner candidate for the 
dark-matter particle [a|- 

In this note, we ask whether resonant coannihilation of 
neutrinos with dark-matter superpartners might be ob- 
servable. We study a particular representative case of 
neutralino dark matter, in which the absolutely stable 
lightest supersymmetric particle (LSP) is Xij a superpo- 
sition of neutral wino, bino, and Higgsino. In circum- 
stances apt for detection, resonant formation of a sneu- 



trino in the reaction vxi 
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shown in Figure ^ might 
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induce absorption lines or direct signatures. A sneutrino 
heavier than the neutralino is implied by the assumption 
that Xi be the LSP. A direct signature of sneutrino forma- 
tion and decay requires that v decay into channels beyond 
the i^Xi formation channel. Informative examples include 
the parameter sets (P and L') presented among the post- 
WMAP benchmarks for the constrained minimal super- 
symmetric standard model in Ref. jTj. (The same two 
models have been considered by Datta and collabora- 
tors [l^l in their recent study of the Lorentz-boosted sit- 
uation of ultrahigh-energy neutralinos scattering on the 
relic neutrino background.) 

We summarize the relevant information in Tabled The 



TABLE I: Superpartner parameters in two MSSM scenar- 
ios [lH|. evaluated using the SPheno p^ code for particle 
properties and the MicrOMEGAs |ll4l code for relic densities, 
as implemented at the Comparison of SUSY spectrum gener- 
ators web site p^ with mt = 172.5 GeV |lal and the default 
values ms = 4.214 GeV and a^Mf) = 0.1172. 
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B{u^ 
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141 GeV 

0.105 

7.7 X 10"^ cm"^ 



< 400 cm"'^' 

265 GeV 
265 GeV 

290 GeV 
299 GeV 
324 MeV 

40.2 nb 
13.0 nb GeV 

0.720 
0.086 
0.194 

277 GeV 
272 GeV 
612 MeV 

52.3 nb 
32.0 nb GeV 

0.342 
0.012 
0.026 
0.620 



185 GeV 

0.111 

6.4 X 10"^ cm"^ 



< 300 cm"^ 

349 GeV 
349 GeV 

422 GeV 

481 GeV 

857 MeV 

20.6 nb 

17.7 nb GeV 
0.434 
0.179 
0.386 

386 GeV 

403 GeV 

2.011 GeV 

14.3 nb 

28.8 nb GeV 

0.153 
0.023 
0.052 
0.772 



first point to notice is that the neutrino energy at res- 
onance is approximately 300 and 500 GeV for the two 
superparticle spectra under study — squarely in the at- 
mospheric neutrino range. Accordingly, there is no need 
to invoke unknown mechanisms to produce the required 
neutrino "beam," as we must for the vv -^ Z^ resonance. 
The modest resonant energies follow from the fact that 
the neutralino mass is more than eleven orders of magni- 
tude larger than the (relic) neutrino mass. In both cases, 
inelastic processes — sneutrino decays that do not return 
to the entrance channel — are prominent. 

There ends the good news, at least for the Universe at 
large. Whereas we expect that stable neutrinos should 
be the most abundant particles in the Universe after the 
photons of the cosmic microwave background, the neu- 
tralino gas is on average very tenuous. The neutralino 
fraction of the Universe (identified with the dark-matter 
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where the numerical value of the critical density is 

Q^ = 1.05 X IQ-^h^ GeV cm~^ (2) 

Consequently the mass density of relic neutralinos is 



g^ = 1.05 X 10 
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• n^h^ GeV 
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= M^M{x'i), (3) 



where A/ij is the mean number density of relic neutralinos 
throughout the Universe. In the two models we consider 
here, N^ ^ 10~® cm~'^, some ten orders of magnitude 
smaller than the relic neutrino density and 31 orders of 
magnitude smaller than the density of electrons in water. 
The peak cross sections for inelastic sneutrino forma- 
tion 
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(4) 
are about an order of magnitude smaller than that [2j for 

Z" formation, a{vv -^ Zfjr^^^^ = Q.ma{w -^ Zy''^ = 
365 nb. The cross sections integrated over the peak 
are some fifty to seventy times smaller than Tz<y{i>D —> 

^)vis'ibic = 912 nb GeV, where Tz = 2.4952 GeV. At the 
resonance peak, the interaction length for neutrinos on 
neutralinos in the present Universe is 



£l, = l/ar'-M^ (5) 

For the two spectra under consideration, we find C^^^^ « 
(1.2,2.4) X lO^^ Mpc for model (r,L'), for ly^^^ interac- 
tions, with similar values for ^'^- interactions. These val- 
ues are roughly 11 orders of magnitude longer than the 
interaction length at resonance for vP —^ Z°, which is 
C\'^^ « 1.2 X 10'' Mpc in the current Universe p|. A dis- 
tance of 10^^ Mpc corresponds to approximately 10^^'^ 
light years, so exceeds the 1.3 x 10^°-year age of the Uni- 
verse by a prodigious factor. If t^Xi ~* ^ coannihilation 
were to be observable, absorption lines would not be the 
signature! We should have to rely on the direct detection 
of few-GeV 7 rays or charged particles from the inelastic 
decay chains that lead back to the :^Xi ground state. 

Our location in the Universe may not be privileged, 
but it is not average. Cold dark matter clusters in galax- 
ies. Most of the analytic expressions proposed in the 
literature for dark-matter halo profiles can be cast in the 
form 
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where go is a characteristic density, a sets the radius of 
the halo, 7 is the asymptotic logarithmic slope at the cen- 
ter, P is the slope as r — > 00, and a controls the detailed 
shape of the profile in the intermediate regions around 
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3.8847 
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3.8351 
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0.0536 
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TABLE 11: Parameters of four radial density profiles of the 
Milky Way dark halo considered in Ref. [l7|l . according to 
the parametrization given by expression ©, plus the mass 
enclosed in a sphere of radius r© . 



Model Q /3 7 a [kpc] go [GeV cm~^] /iQ [10®^ GeV] " 



ISO 2 2 4.0 

BE [18j 13 0.3 10.2 

NEW [19j 1 3 1 16.7 

M99 [201 1.5 3 1.5 29.5 



a. A useful catalogue of dark-matter halo profiles for the 
Milky Way galaxy is given in Ref. [i3|- The parameters 
of four profiles examined there are collected in Table E] 
In each case, the density parameter Qq is adjusted to re- 
produce the presumed dark-matter density (neglecting 
local influences) at our distance (r© = 8. 5 kp c) from the 
galactic center, g^(r0) = 0.3 GeV cm~'^ [21|, which rep- 
resents five orders of magnitude enhancement over the 
relic density in the Universe at large. The resulting neu- 
tralino number densities, Af^{rQ) = g^{rQ)/M{xi), are 
given in TableHJ the density profiles themselves are shown 
in Figure |21 along with the dark mass enclosed within a 
sphere of specified galactocentric radius. 

What little we know about the distribution of dark 
matter within the solar system has been deduced by in- 
vestigation possible dark-matter perturbations on the or- 
bits of planets and asteroids. In a recent contribution, 
Khriplovich and Pitjeva Q^ have deduced bounds on the 
amount of dark matter that lie within the orbits of Mer- 
cury, Earth, and Mars. On the assumption that the dark- 
matter density is constant within the spherical volumes 
delimited by these three planetary orbits, the Earth and 
Mars orbits constrain 



P|<10 
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g cm-3 = 5.6 X 10^ GeV cm^^^ (7) 



which would imply upper bounds on the number density 
of neutralinos at Earth, J\f^ < 350 cm^"^. 

What are the prospects for observing the coannihi- 
lation of neutrinos on the relic neutralino background? 
Consider first the interaction of neutrinos created in 
the collisions of cosmic rays with Earth's atmosphere. 
According to the calculations of the atmospheric neu- 
trino flux reviewed by Gaisser & Honda |23|, the muon- 
neutrino flux at E^ « 300 GeV is 
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The quantity Fj; cr-.'jnoi provides a rough measure of the 



coannihilation cross section integrated over the resonance 
peak. Consequently, the product 
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j^ peak 



dN, 
dE^ 



(9) 




10-' 10° 10' 

r (kpc) 
FIG. 2: Top panel: Dark-matter density profiles given by 
the form @ with parameters specified in Table ^ The 
solid (blue) curve shows the canonical Navarro-Frenk-White 
universal profile fl9| . The dot-dashed (purple) curve repre- 
sents the cuspier profile of Moore and collaborators !i20| . Also 
shown are the Binney-Evans profile 18] (dotted green line) 
and the nonsingular profile corresponding to an isothermal 
sphere (dashed red line). Bottom panel: The corresponding 
enclosed-mass profiles. 



characterizes the v coannihilation rate of a single neu- 
tralino in the rain of atmospheric neutrinos. 

For model F, we compute the coannihilation rate 

TZ}^ «4.3xl0~''°s^^. If we consider the neutralino 
density for the Universe at large, then the rate per unit 
target volume is 



nP^M<{R 



(I') 



3.3 X 10"''^ cm~3 s"^ 



10 



-« cm-3 y-\ 



(10) 



which is laughably small, and smaller still for model L'. 
The rates (jTO)) are multiplied by (2.9 x 10^,5.2 x 10^") 
if we consider the mean density at Earth's distance from 
the galactic center or the orbital-dynamics upper bound 
(|7|) on the density of dark matter in the solar system. 

To pursue this line to a logical conclusion, we ask how 
many interactions of atmospheric neutrinos with relic 
neutralinos occur in the atmosphere per unit time. Tak- 
ing the Earth to be a sphere with radius R^ = 6 371 km, 
we compute the volume of a 10-km-thick shell above the 
Earth's surface to be T4tm = 5.1 x 10^** cm'^. The num- 
ber of sneutrino events induced in this volume is therefore 
^ 3 X 10~^ y~^, taking the largest conceivable neutralino 
number density. Such an infinitesimal rate renders moot 
a discussion of signatures and detection efficiencies. Set- 
ting aside questions of detection, the total number of 



events within the Earth's volume, I/q = 1.08 x 10^"^ cm'^, 
is no more than 6 x 10~^ y~^. 

Extraterrestrial sources may interact with relic neutri- 
nos over a larger volume. A representative estimate |24| 
of the diffuse neutrino flux from active galactic nuclei 
suggests that, at E^ « 300 GeV, the flux of cosmic neu- 
trinos arriving from all directions is 



dE,, 



6.3 X 10"" cm-2 g-i Qev^i 



(11) 



approximately 2 x 10"^ x the vertical atmospheric neu- 
trino flux we have just considered. The upper bound 
on the amount of dark matter in the sphere defined by 
Earth's orbit, //(I au) < 7.85 x 10'*'' GeV ^, could 
therefore be the site of up to (144,130) z^Xi ~* ^ inelastic 
interactions per year for model (r,L'). 

Let us compute the contribution of each event, at its 
own location, to the signal recorded by a detector of unit 
area at Earth. It will be sufficient to assume that the 
detector records signals from all relevant directions with 
perfect efficiency. We define the vector from the Sun 
to the detector as s = (0,0, s), and denote the location 
of the target as f = r(sin6'cos0, sin^sin^, cos6'). Then 
the vector that points from the target to the detector is 



d = s — f, so that d 



2rscos9. The decay 



products of the sneutrino produced at r are distributed 
isotropically. Accordingly, the fraction incident on a de- 
tector of area A will be / = A/ATrd"^, the fraction of the 
solid angle that the detector subtends. 

The effective number of targets seen by a detector of 
unit area at s is given by 



n-off(s) 



1 



47rAf-f 



'3„ 6'x('^) 
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r^ drQj^{r) 



p2 _|_ g2 _ 2rs cos 6 



For the special case of a constant density g^ = 5.6 x 
10"* GeV cm-3, and s = 1 au = 1.496 x lO^^ ^m, the 
effective number of targets is 



neff(s) 



2M^ 



3 x 10 



15 cm-2. 



(13) 



Using the cosmic- neutrino flux (|ll|l . we estimate that a 
detector in the vicinity of Earth would be sensitive to 
7.7 X 10^26 events cm-2 y-i^ With such a small num- 
ber of events, there is no hope of detecting the few-GeV 
gamma rays that would be created in the cascade back 
to X?. 

The halo of our galaxy contains a great quantity of 
dark matter: for the NEW proflle, 4.4 x 10^^ GeV hes 
within the galactocentric radius of our solar system. The 
number of neutralino targets with which neutrinos might 
coannihilate to form sneutrinos is thus (3.1,2.4) x 10^^ 
for model (T, L'). Focusing again on model F and taking 
the estimate (|ll|l for the cosmic-neutrino flux, we find 
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FIG. 3: Effective number of neutralino targets per unit of 
detector area viewed from a location at radius s from the 
galactic center, according to four dark-matter density profiles 
for the Milky Way galaxy. Our location corresponds to s = 
8.5 kpc. 



that the sneutrino formation rate throughout the galaxy 
is 2.6 X 10*^ s"* = 8.1 X 102^ y"*. These rates are prodi- 
gious in absolute terms, but represent an insignificant 
disturbance to the galaxy's neutralino population. 

Can we hope to observe the sneutrino formation that 
might be bubbling away in our neighborhood? in this 
case, we define s = (0, 0, s) to be the vector from the 
galactic center to the detector, and use (I12|) to compute 
the effective number of targets seen from Earth. For the 
NFW profile, we find ncs{rQ) = 4.9 x 10"^^ cm-2. p^j, 
a detector of any plausible area, only a tiny fraction of 
the 3.1 X 10^5 neutralinos in the halo produce a visible 
signal. Indeed, using the neutrino fiux l|ll|l . we expect 
that a detector in the vicinity of Earth would be sensitive 
to 1.3 X 10"2i events cm-2 y-i 

It is not obvious that Earth's location should happen 
to be optimal for viewing coannihilation events in the 
galaxy. In Figure |31 we show how the effective number 
of targets viewed depends on s, the galactocentric radius 
of the observation point. For the smooth ISO and BE 
profiles, our solar system lies near the optimal distance, 
while for the NFW profile the effective number of tar- 
gets viewed is insensitive to the detector position. The 
most singular profile we consider, M99, favors an obser- 
vatory close to the galactic center — not that we could 
contemplate one — but even in this case, the sensitivity is 
enhanced by less than two orders of magnitude, which is 
far too little to enable detection. Singular profiles have 
more effect on the rates for neutralino-neutralino anni- 
hilations, which are proportional to the square of the 



neutralino density. 

If ncutralinos account for much of the cold dark mat- 
ter of the Universe, then it is possible that neutrino- 
neutralino coannihilation into sneutrinos is happening all 
around us, but at a rate that will forever elude detection. 
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